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^ ■ Abstract 

o : 

^ ; We study the fate of U(l) strings embedded in a non- Abelian gauge theory 

with the hierarchical pattern of the symmetry breaking: G — — > U(l) — nothing, 
V ^> v. While in the low-energy limit the Abrikosov-Nielsen-Olesen string (flux 
tube) is perfectly stable, being considered in the full theory it is metastable. We 
consider the simplest example: the magnetic flux tubes in the SU(2) gauge theory 
with adjoint and fundamental scalars. First, the adjoint scalar develops a vacuum 
expectation value V breaking SU(2) down to U(l). Then, at a much lower scale, 
the fundamental scalar (quark) develops a vacuum expectation value v creating 
the Abrikosov-Nielsen-Olesen string. (We also consider an alternative scenario in 
which the second breaking, U(l) — — > nothing, is due to an adjoint field.) We 
suggest an illustrative ansatz describing an "unwinding" in SU(2) of the winding 
inherent to the Abrikosov-Nielsen-Olesen strings in U(l). This ansatz determines 
an effective 2D theory for the unstable mode on the string world-sheet. We calculate 
the decay rate (per unit length of the string) in this ansatz and then derive a general 
formula. The decay rate is exponentially suppressed. The suppressing exponent is 
proportional to the ratio of the monopole mass squared to the string tension, which 
is quite natural in view of the string breaking through the monopole-antimonopole 
pair production. We compare our result with the one given by Schwinger's formula 
dualized for describing the monopole-antimonopole pair production in the magnetic 
field. 



1 Introduction 



Formation of the chromoelectric flux tubes in non-Abelian gauge theories is being 
discussed as a mechanism of color confinement, at a qualitative level, since the early 
days of quantum chromodynamics (QCD). No quantitative first-principle description 
of the phenomenon was ever constructed in QCD — strong coupling regime inherent 
to this theory precluded all efforts in this direction. 

The revival of interest to non-Abelian gauge theories occurred in the mid-1990's. 
The idea of using supersymmetry as a tool allowing one to deal, to a certain extent, 
with strong coupling regime was revived. This development culminated in the work 
of Seiberg and Witten |TJ, who considered M = 2 Yang-Mills theory slightly 
perturbed by a (small) mass term of the adjoint matter field. This perturbation 
breaks M = 2 down to M = 1. Supersymmetry proved to be sufficiently powerful 
to allow Seiberg and Witten to constructively demonstrate the existence of the dual 
Meissner effect — the monopole condensation accompanied by the formation of the 
chromoelectric flux tubes. Technically, the Seiberg- Witten theory has two distinct 
scales (this was crucial for their construction): the scale of strong interaction A, and 
a much smaller scale regulated by a small adjoint mass term. Below A the original 
SU(2) Yang-Mills theory reduces to an Abelian (dual) quantum electrodynamics 
(QED). Correspondingly, the flux tubes (strings) of the Seiberg- Witten theory are 
the conventional U(l) Abrikosov-Nielsen-Olesen (ANO) strings |3], |], ||. In the 
SU(A^) case, the low-energy effective theory is that of U(l) Ar_1 ; the flux tubes one 
deals with present a straightforward generalization of the ANO strings. 

In the Yang-Mills theory without quarks, the underlying gauge group has a non- 
trivial 7Ti homotopy. On topological grounds one then expects Zn vortices to appear 
H B H; H> 13 01 • However, in fact, in the Seiberg- Witten construction, in which 
at low energies the gauge group is broken down to U^l)^ -1 by the adjoint matter 
condensation, it is the U(l) Ar_1 strings that appear [12j near the monopole/dyon 
vacua. They form N — 1 infinite towers of the ANO flux tubes for each of U(l) factors 
and give rise to confinement of quarks. It is clear that only some of these strings 
(those which correspond to Zn strings of the microscopic non-Abelian theory) are 
stable, all others must be unstable (metastable) [13[] beyond the extreme low-energy 
limit. 

A slightly different scenario takes place in the Seiberg- Witten theory with fun- 
damental matter (quarks). In this case the underlying gauge group SU(A r ) has a 
trivial homotopy group, 7Tl(SU(./V)) = 0, and does not admit flux tubes. However, 
near the charge vacua this theory has an Abelian low-energy description too, which 
ensures the presence of the ANO flux tubes in the extreme low-energy limit. These 



ANO strings give rise to confinement of monopoles JTJ, [15|, [16], \Fjj. It is perfectly 
clear that all these ANO strings must be metastable in the microscopic non-Abelian 
theory. 

Thus, we see that the following general question presents a considerable interest. 
Assume that one considers an underlying non-Abelian gauge theory with a trivial (or 
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"almost trivial") homotopy group 7Ti. (For definiteness, one may choose SU(2).) As- 
sume that this theory experiences a two-stage spontaneous symmetry breaking: first, 
at a high scale V, the non-Abelian group is broken down to an Abelian subgroup (let 
us say, U(l), for definiteness), and then this U(l), in turn, is spontaneously broken 
at a much lower scale v. The underlying non-Abelian theory will be referred to as 
"microscopic," while the low-energy Abelian theory will be called "macroscopic." In 
the low-energy limit one can forget about the microscopic theory and consider QED 
with a charged matter field which develops a vacuum expectation value (VEV). 
Correspondingly, stable ANO strings do exist in the macroscopic theory. In the 
microscopic theory they are metastable, rather than stable, however, because of the 
triviality of tti of the original gauge group. A long ANO string can (and will) break 
due to the monopole-antimonopole pair creation. The task is to find the probability 
(per unit length per unit time) of the string breaking. 

This question can be quantitatively addressed at weak coupling. In this paper we 
consider a simple (non-supersymmetric) model which closely follows the pattern of 
the (supersymmetric) Seiberg-Witten theory. We start from an SU(2) gauge model 
with scalar fields — one in the adjoint and another in the fundamental represen- 
tation — with a certain interaction between them. The field in the fundamental 
representation will be referred to as the "quark field." The interaction of scalars is 
arranged in such a way that the adjoint scalar develops a large vacuum expectation 
value, 

^>A, (1.1) 

where A is the dynamical scale of the SU(2) theory. This VEV of the adjoint field 
breaks the SU(2) gauge group down to U(l) and ensures that the theory at hand 
is weakly coupled. At this stage the 't Hooft-Polyakov monopoles emerge |1| |I9| . 
Their mass is very heavy, Mm ~ V/g. Below scale V one is left with QED. The 
(charged) quark field develops a small VEV v, 

v<^V. (1.2) 

Then the standard ANO flux tubes emerge. We study their decay in the qua- 
siclassical approximation. To this end we consider dynamics of an unstable mode 
associated with the possibility of "unwinding" the ANO string winding on the SU(2) 
group manifold. This is an under-barrier process, with the corresponding action be- 
ing very large in the limit v -C V. The physical interpretation of this tunneling 
process is the monopole-antimonopole pair creation accompanied by annihilation of 
a segment of the string. 

We present an analytic ansatz which explicitly "unwinds" the string. The string 
decay rate is found in the framework of this ansatz. Our task — constructing fully 
analytic ansatz — is admittedly illustrative. Therefore, we limit ourselves to a 
minimal number of profile functions. The ansatz obtained in this way is suitable 
for a qualitative understanding of the phenomenon. It is too restrictive to describe 
the production of the monopoles per se; rather, it describes the production of a 
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pair of highly excited states with the monopole (antimonopole) quantum numbers. 
We argue, however, that an effective description of the string decay suggested by 
this ansatz is more general and is applicable for realistic ansdtze provided that (i) 
the condition ( |1.2|) is met, and (ii) the string tension and the monopole mass are 
treated as given parameters. In terms of these parameters the general result for the 
metastable string decay rate (the probability per unit time per unit length of the 
string) is 

2 f 7T Mh ) . . 

Tbreaking ~ V exp ^ - — \ , (1.3) 

t -L ANO J 

where Mm is the monopole mass and Tano is the string tension. It is worth remind- 
ing that M 2 M ~ V 2 /g 2 while T AN o ~ v 2 . 

As was expected, Tweaking turns out to be exponentially suppressed. This is nat- 
ural in light of the tunneling interpretation of the string breaking process — that 
the string is broken into pieces by the monopole- ant i- monopole pair production. In 
fact, our result is quite similar to Schwinger's formula pOj dualized to the mag- 



netic charge production in the constant magnetic field [pi]] . The only difference is 
the replacement of Bg~ 2 in the dualized Schwinger formula by Tano- This is not 
surprising: the macroscopic theories of these two phenomena are similar. 

The result presented in Eq. ( |1.3| ) is not new. It was obtained long ago in Refs. 



23| . In these papers the decay of metastable ANO strings arising in theories with 
a hierarchical pattern of symmetry breaking was treated in the framework of an 
effective theory, by calculating the action of the bubble formed by the monopole 
world line on the string world sheet. 

In our present paper the logic is different. We "forget" for a while about 
monopoles per se, and use an "unwinding" ansatz to derive a quasi-modulus the- 
ory for the unstable string mode on the string world sheet. We then consider the 
bubble creation in this quasi-modulus theory and rederive Eq. ( |1.3|) . After this is 
done, we interpret ( |1.3| ) as the probability of the monopole-antimonopole pair pro- 
duction. From this standpoint our results can be viewed as a demonstration of the 
assertion that the "unwinding" of a metastable string goes via the production of 
the monopole-antimonopole pairs. We thus provide a necessary background for the 
effective approach formulated in the pioneering works [TZ^, [23|] . 

If the small- VEV field of the quark type (i.e. in the fundamental representation) 
is replaced by a small- VEV field in the adjoint (such that the VEV's of two adjoint 
fields present in this model are misaligned), then one arrives at a weak coupling 
model of Z 2 strings. The gauge symmetry of the microscopic theory is now SU(2) / 
Zi- Since 7Ti(SU(2)/Z2) = Z% , the minimal magnetic flux string is absolutely 
topologically stable. Higher-flux strings are stable only in the low-energy U(l) limit, 
and can decay through the monopole-antimonopole pair production. We construct 
an "unwinding" ansatz in this case too, and argue that the decay rate is given by 
the same expression (|1.3|). 

The paper is organized as follows. In Sect, ^we formulate our model and explain, 
in concrete terms, what needs to be done in order to calculate the string decay 
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probability. In Sect. f|we present an ansatz for the gauge and scalar fields describing 
the under-barrier transition ( "unwinding" ) under consideration. The decaying string 
ansatz is parametrized by three profile functions. We identify an unstable mode 9 
in which tunneling occurs. For the extreme type II and type I strings one needs to 
know only the asymptotic behavior of the above profile functions. This simplifies 
the problem immensely. 

Section [5] is devoted to the extreme type II string which arises in the theory 
with the quark mass much larger then the photon mass (m q 3> m 7 ). In this limit, 
the stable ANO solution was obtained analytically by Abrikosov long ago ||. We 
generalize this solution to cover the decaying string. For the unstable mode 9 we 
derive an effective two-dimensional field theory on the string world-sheet. The string 
breaking in the microscopic theory corresponds to the false vacuum decay in the 
effective world-sheet field theory for 9{t, z). We then apply well-developed methods 
24|, 2~5 for calculating the false vacuum decay rate through bubble formation. In 



this way we find the probability of the string breaking for the extreme type II strings. 

In Sect. H we consider the opposite case of the extreme type I string {m q -C m 7 ). 
The analytic solution in this case was obtained quite recently in Ref . |0| . Although 
the string tension for type I is given by an expression significantly different from that 
for type II, the string decay rate, being expressed in terms of the string tension, is 
determined by the same formula as in the type II case. 

In Sect. 7 we turn to the theory with two adjoint scalars and work out the decay 
rate of the string with winding number n = 2. Note, that the string with minimal 
winding number n = 1 (^-string) is stable in this model. Section [8] extends the 
proof of Eq. ( |1.3| ) to non-extreme strings. 

In Sect. P we compare our result for the string decay rate with the one given by 
Schwinger's formula [20] which might be used for the evaluation of the probability 
of the monopole-antimonopole pair production in the external magnetic field. The 
comparison can be performed only qualitatively. 

The reason is that Schwinger's formula and Eq. ( \1.3\) refers to different phases of 
the theory. Schwinger's formula deals with the monopole production in the external 



magnetic field in the Coulomb phase, while Eq. (|1.3|) describes breaking of ANO 
string by monopole-anti-monopole pair in the phase in which monopoles are confined. 
Still, the comparison exhibits a qualitative agreement with result ( |1.3|) as far as the 
powers of the mass scales in the exponent are concerned. 



Finally, Sect. 10 summarizes our results and conclusions and outlines problems 



for future investigation. Appendix contains some details of an "improved" unwinding 
ansatz. 



4 



2 The model and formulation of the problem 



In the bulk of the paper we will consider SU(2) gauge theory with the action 

(2.r 




(a = 1,2,3) is a real scalar field in the adjoint, while (k = 1, 2) is 
a complex scalar field in the fundamental (sometimes, we will refer to it as to the 
"quark" filed). Finally, g is the gauge coupling, and V(q, </>) is a scalar self-interaction 
potential. Throughout the paper we will deal with the adjoint fields both, in the 
matrix and vector notations, say 

The covariant derivatives and V M act in the adjoint and fundamental repre- 
sentations, respectively. The simplest form of the potential V(q, (p) that will serve 
our purpose is 



v( q , 0) = a (| g | 2 - v 2 ) 2 + a (rr -v 2 ) 2 + 1 1 (V - ^) 



Q 



(2.2) 



where v and V are parameters of dimension of mass and A , A and 7 are dimensionless 
coupling constants. In this work we limit ourselves to the case of weak couplings, 
when all four coupling constants g 2 , A, A and 7, are small. We also assume that 
V ^> A, where A is the scale parameter of the SU(2) gauge theory. Then the 
quasiclassical treatment applies. Since our goal is a non-perturbative string decay, 
we will ignore perturbative quantum corrections altogether. 

To arrange the double-scale (hierarchical) pattern of the symmetry breaking 
mentioned in Sect. |l] we must ensure a hierarchy of the vacuum expectation values 
(VEV's). Namely, the breaking SU(2)— >U(1) occurs at a high scale, while U(l)— > 
nothing at a much lower scale, 

v<&V. (2.3) 

At the first stage the adjoint field develops a VEV which can be always aligned 
along the third axis in the isospace, 

((j) a ) = 5 a3 V. (2.4) 

This breaks the gauge SU(2) group down to U(l) and gives masses to the 
bosons, and to one real adjoint scalar </> 3 , 

m w ± = gV , m adj = m a = 2\f2~\V , (2.5) 

while two other adjoint scalars ((f) 1 and <fi 2 ) are "eaten" up by the Higgs mechanism. 
Note that simultaneously the second component of the quark field, 92, acquires a 
large mass, 

M q2 = ^V, (2.6) 
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due to the last term in the potential (p.2|). 

Below the scales (|2.5|), ( |2.6|) the effective low-energy theory reduces to QED: the 
U(l) gauge field interacting with one complex scalar quark q 1 . The action is 



5qed= / d'x l± F^F^ + V, qi - + \{\ qi \ 2 -vJ\ , (2.7) 



2 . /, , 2 2 X 2 



where 



Furthermore, at this second stage the charged field q± develops a VEV, and the 
U(l) theory finds itself in the Higgs phase, 

(qt) = v, (while (q 2 ) = 0). (2.8) 

At this stage the gauge group is completely broken. The breaking of U(l) gives a 
mass to the photon field A 3 ^, namely, 



1 / N 

m 1 = ^gv, (2.9) 



while the mass of the light component of the quark field q± is 



m, 



2\f\v. (2.10) 



In what follows we will essentially forget about the heavy component of the quark 
field q2, it will be irrelevant for our consideration 0. Only qi is relevant. Correspond- 
ingly, in the bulk of the paper we will drop the subscript 1 in mentioning the quark 
field; by definition, m q = m qi = 2yXv. 

The theory ( |2.7| ) is an Abelian Higgs model which admits the standard Abrikosov- 
Nielsen-Olesen (ANO) strings J3], |J. Let us briefly review their basic features. For 
generic values of A in Eq. ( |2.7[ ) the quark mass m qi (the inverse correlation length) 
and the photon mass m 7 (the inverse penetration depth) are distinct. Their ratio 
is an important parameter in the theory of superconductivity, characterizing the 
type of superconductor. Namely, for m qi < m 7 one deals with the type I supercon- 
ductor in which two strings at large separations attract each other. On the other 
hand, for m qi > m 7 the superconductor is of type II, in which two strings at large 
separations repel each other. This behavior is related to the fact that the scalar 
field generates attraction between two vortices, while the electromagnetic field gen- 
erates repulsion. The boundary separating superconductors of the I and II types 
corresponds to m qi = m 7 , i.e. to a special value of the quartic coupling A, namely, 

A = — . (2.11) 



lr The only place where qi surfaces again, implicitly, is in Eq. (4.1) at 9 ^ 0. There is no menace 
of confusion, however. 
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In this case the vortices do not interact. 

It is well known that the point ( 2.11|) represents, in fact, the Bogomolny-Prasad- 
Sommerfield (BPS) limit. At m qi = m 7 the ANO string satisfies first order differen- 
tial equations and saturate the Bogomolny bound [Q. In supersymmetric theories 
the Bogomolny bound for the BPS strings coincides with the value of the central 
charge of the SUSY algebra [^6|, ^7], ^8|. In particular, the BPS strings arise in 
the Seiberg-Witten theory near the monopole/charge vacua at small values of the 
adjoint mass perturbation [TJ], |29], |50| . 

As was mentioned in Sect, [l], we will also consider, in brief, a model with two 
adjoint matter fields, 

S = J d*x^F^ a + \{D,r) 2 + \{D, X a Y + U(x^)\ , (2.12) 

where a and x a (o = 1, 2, 3) are real scalar fields, and the potential U(x, 0) can be 
chosen, for instance, as follows: 

U(xA) = a(r<p a -V 2 ) 2 + ^ (x a X a ~ 2v 2 ) 2 + 1 (<p a X a ) 2 (2-13) 

The condensation of the field, Eq. ( |2T4]) , breaks the gauge symmetry SU(2)/Z 2 — > 
U(l), and makes W bosons, 3 and x 3 heavy. What remains in the low-energy limit? 
The low-energy is very similar to ( |2.7| ), namely, 



S QED , = J ^ |_L F 3^^3 + + p ( x+x - _ ) (2M) 



where 



1 

X' 



± " 4? x*±ix 2 



- V2 

= (d,±iAl) X ± . (2.15) 

At the second stage the remaining U(l) is completely broken by the condensate of 
the x field, 

\(x)\=v. (2.16) 
This gives mass to the "photon" field A 3 , and the ANO strings enter the game. 



3 Abrikosov-Nielsen-Olesen String 



In the model ( |2.7f ) the classical field equations for the ANO string with the unit 
winding number are solved in the standard ansatz, 



qi{x) = q{r) e %OL 



A 6 



A\{x) = 2 £ii ^ [1 - f(r)} . (3.1) 
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Here 

r = ,/e a 

is the distance from the vortex center while a is the polar angle in the transverse 
to the vortex axis (l,2)-plane (the subscripts i,j = 1,2 denote coordinates in this 
plane, x and y, see Fig. |1]). Moreover, q(r) and /(r) are profile functions. Note, 
that did = —€ij Xj/r 2 . 

t Z 




Figure 1: Geometry of the string. 

The profile functions q and / in Eq. ( |3.1| ) are real and satisfy the second order 
differential equations 



q + -q - —f q -m — = 

r r Zv 

f" ~ -f'--iq 2 f = (3.2) 



for generic values of A (the prime stands here for the derivative with respect to r) 
plus the boundary conditions 
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q(0) = 0, /(0) = 1, 

q(oo)=v, /(oo) = 0, (3.3) 

which ensure that the scalar field reaches its VEV (q\ = v ) at infinity and the vortex 
at hand carries one unit of magnetic flux. 

The expression for the string tension (energy per unit length) for the ANO string 
in terms of the profile functions ( p.lQ has the form 

T ANO = 2njrdr j A £ + g' 2 + £f + A(g 2 - , 2 ) 2 } . (3.4) 

For generic values of the ratio m 9 /m 7 only a numerical solution of Eqs. ( |3.2| ) 
is possible. However, in the extreme type II case (m g ^> m 7 ) and extreme type I 
case {m q <C m 7 ) analytical solutions can be readily found ||, |i~5j . We will review 
these solutions in Sects. || and |6], respectively. As was explained in Sect. |], in the 
full SU(2) theory the ANO string can decay. We will use the solutions discussed in 
Sects. H and |6] in order to analytically calculate the decay rate of the ANO string. 

The magnetic field flux for the string (|3.1|) is 

- J B 3 dxdy = ^ j A 3 i dx i = 27i . (3.5) 

In the model with two adjoints, Eq. ( |2.14j ), the ANO string solution has the form 

X' = h /2 (r)e- ia , /i(0) 1/2 = 0, h 1/2 (oo) = v , 
Al = 0, 

A\{x) = e i3 ^ [1 - f(r)] , /(0) = 1, /(oo)=0. (3.6) 

In this case 

J B 3 dx dy = j A\dxi = 2tt . (3.7) 

The magnetic flux is twice smaller. That's the reason why the elementary string in 
the model (|2.14j ) cannot be broken by the monopole pair production. However, the 
double-winding string 

x - = h(r)e-* 2a , 
A 3 = 0, 

A]{x) = 2e ij ^[l-f(r)\ (3.8) 

can and will be broken. 
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4 Decaying strings: a " primitive" unwinding 
ansatz 

It is clear that the ANO strings ( |3.1| ) are topologically stable only at low energies 
when the SU(2) theory Q2.1]) reduces to a "macroscopic" theory, QED, see Eq. (|2.7| ). 
In the full "microscopic" theory Q2.ip they should be metastable because the SU(2) 
gauge group does not admit flux tubes, 7Tl (577(2)) = 0. To visualize the decay 
possibility, note that the winding in ( |3.1| ) runs along the "equator" of the SU(2) 
group space (which is S3) and, therefore, can be shrunk to zero by contracting the 
loop towards the south or north poles (Fig. g). 

t 3 




SU(2) group space 



Figure 2: Unwinding the ANO ansatz. 

It is not difficult to devise an ansatz encoding the possibility of unwinding the 
field configuration ( |3.1| ) through the loop shrinkage in the SU(2) group space. The 
ansatz which does the job and eventually will allow us to calculate the ANO string 
decay rate is parametrized by an angle parameter 9, 

q k (x) = U e k q e {r) , {e x = 1, e 2 = 0) , 

A = 0, A 3 = 0, 
Aj{x) = Wd j U- 1 [\-f 9 (r)\, J = 1,2, 

cf) = VU^-U- 1 + (4.1) 
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where 

U = e~ iaT3 cos6 + i n sin 9 . (4.2) 

(Eventually, upon quantization, 9 will become a slowly varying function of z and t, 
afield 6(t,z).) 

The gauge and quark fields in ( |4.1| ) are parametrized by profile functions /e(r) 
and ge(r) depending on the parameter 9. They satisfy the same boundary conditions 



qe(0) = 



oo 



fe(0) 
fe(oo) -- 







(4.3) 



as in the U(l) case, see Eq. ( p.3[ ). The boundary conditions at zero are chosen to 
ensure the absence of singularities of our ansatz at r = 0. The magnetic flux of the 
string ( [4.i|) is 

2ttcos 2 #. (4.4) 

It equals to 27r at 9 = and goes to zero at 9 = |. 

The term A<p in the last line of Eq. (4.1) is needed to make sure that there is no 
singularity at r —>■ 0. For axially symmetric string the function A<p can be chosen 
in the form 



Acf) = <p g (r) 



— sin ot cos a 

2 2 



(4.5) 



where we assume that the component of A0 along t% is zero, while </?e(r) is an extra 
profile function, which depends on 9 as a parameter. The a = 1, 2 components of 
A(j) cannot be put to zero. To see this substitute Eqs. (|4.2|) and ( |4.5|) into the last 
line in Eq. (|4. 1|) . Then one gets 



73 



V cos 29 



— sin a cos a 

2 2 



[Vsm29-ve(r)} . 



(4.6) 



From this expression it is clear that has no singularity at r = provided that 

<p g (0) = V sin20 . (4.7) 
The boundary condition for (fe(r) at infinity should be chosen as follows: 

W(oo) = 0. (4.8) 
Both boundary conditions are consistent with the initial condition 



M r )\ 



0=0 







(4.9) 



to be imposed. 

For future reference it is convenient to present the very same ansatz ( |4.1| ) in the 
singular gauge, 

Qk(x) = e k q e (r), 

A(x) = -iidiU'^Ufeir), 

(j) = V— +U- 1 (A(f))U. (4.10) 
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Our ansatz, Eqs. ( |4.1|) or ( 4.10 ), plus ( |4.2|) , smoothly interpolates between the 
ANO-type winding along the equator at # = 0, and constant matrix with no winding 
at 9 = it/2. In other words, we start from the ANO string at 9 = and arrive at 
empty vacuum at 9 = it/2. Indeed, at 9 = the adjoint field <ft is aligned, <p a = VS 3a , 



while U = e~ iaT3 . At 9 = n/2 again <p a = VS 3a , and U — iri, implying Aj = 0. 
(Equivalently, one could have chosen to end up with U = ir 2 at 9 = |.) Thus, 
we managed to unwind the ANO winding through extra dimension of the vacuum 
manifold in the SU(2) theory which was not there in QED. 

We pause here to make an additional comment regarding our ansatz (4J.). At 
large r, when q 9 — > v and fg — > and <pg — > 0, our field configuration presents a 
gauge-transformed "plain vacuum". This ensures that at every given 9 the energy 
functional converges at large r. The convergence of the energy functional at small 
r is guaranteed by the boundary conditions qe(0) = , fe(0) = 1 and (|4~7| ). 

The tension of the string (or the field configuration in which it evolves at 9 ^ 0) 
is a functional of three functions fe(r) , qo{r) and (pe(r). It is easy to get this 
functional by substituting the ansatz ( |4.1|) in the action ( |2.1|) . Restricting ourselves 
to (l,2)-plane we obtain after some algebra 

T(0) = 2nJ rdr jl S cos 2 * + ^ + S± g 2 cog 2 ^ 



1 



,'2 



+ 



2r 2 



(^ (cos 29 - 2f B cos 2 0) + Vf e sin 2^ 



+ KTe-v) +X tpi-2Vtp e sm29 



2 

+ 4 



2„2 



(4.11) 



Needless to say that at 9 = the string tension T coincides with that for the 
ANO string, see (|3.4| ), while at 9 = -n/2 it goes to zero, as was expected. 

Now we have to minimize the string functional ( |4.11| ) with respect to three profile 
functions fe(r) , qe{r) and (fe(r) at fixed 9. This procedure would give us a solution 
for the profile functions. Finding the full solution is a rather complicated task 
requiring numerical computations which go beyond the scope of the present paper. 
However, we will be able to get sufficient insight in order to derive the general 
formula ( |1.3|) by purely analytical means. 



Why we call the ansatz (|4. 1| ) primitive? It has only one profile function fo(r) 
for all three gauge components. This is okay at 9 = when the W boson degrees 
of freedom are not excited. At 9 ^ the universality of fe{r) forces one and the 
same spread in the perpendicular plane of the W boson and photon components of 
the solution. As a result, unwinding of the ANO strings in Eq. (|4.1|) proceeds via 
the production of highly excited "monopoles." This shortcoming can be eliminated 
through introduction of extra profile functions, see Sect. |8] and Appendix. 
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5 Extreme type II string 



As we have already mentioned, even for the stable ANO string the analytic solu- 
tion for generic values of the ratio m q /m 1 is absent. On the other hand, it is not 
difficult to find the solution that would be valid in the logarithmic approximation, 
lnm 9 /m 7 ^> 1. Therefore, in this section we consider the extreme type II supercon- 
ductor for which we will be able to obtain an analytic solution of the problem of the 
decaying string. 

In fact, we impose the following relation between the adjoint scalar, quark and 
photon masses: 

m a > m q > m 7 . (5.1) 

In addition, it is convenient to limit ourselves to the case m^y ~ m a . At first, 
we review Abrikosov's solution for the ANO string ||, then consider the decaying 
string and, finally, work out the effective action for the unstable mode on the string 
world-sheet and calculate the string decay rate. 



5.1 Type II ANO string 

Assume that m q ^> m 7 in the Abelian Higgs model ( |2.7| ). Then the ANO string 



looks as follows. The quark field q(r) varies from zero to its vacuum value v inside 
a small core of radius of order of l/m q , whereas the electromagnetic field is spread 
over a much larger domain, of order of l/m 7 . In the latter domain the quark field 
is already very close to its VEV. The solution of the second equation in (|3.2|) for 
the gauge profile function / is / ~ r 2 in this domain (where the last term in this 
equation can be ignored). Being properly normalized, the solution has the following 
asymptotics: 

!1 — emir 2 at r <C m~ l , 
(5-2) 
at r > m 7 , 

where c is a constant, c ~ 1. The approach to zero at r ^> m^ 1 is exponentially 
fast. 

The leading (and the only) logarithmic contribution to the string tension comes 
from the third term in the expression ( |3.4| ) for the ANO string tension. In the domain 
m" 1 C r <C rn~ l , with the logarithmic accuracy, one can substitute in Eq. (3.4|) 
q = v and retain only the leading term in the gauge profile function, f — 1 + ■ ■ ■. In 
this way one gets 

T%o = 2™ 2 In -i . (5.3) 

m 7 

Domains other than m^ 1 < r <C m^ 1 yield a non-logarithmic contribution. The 



same is valid with regards to the first, second and fourth terms in Eq. ( |3.4| ), as 
well as deviations from q = v and / = 1. All these effects give corrections to ( |5.3|) 
suppressed by powers of 1/ In (m q /m 7 ). 
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This result was obtained by Abrikosov in 1957 ||. Note that the leading log- 
arithmic contribution to the string tension ( |5.3|) is totally insensitive to details of 
the profile functions q and /. It "feels" only the boundary values q(oo) = v and 
/(0) = 1. 

5.2 Decaying type II string 

Now we use the same method of separating distinct physical scales to describe the 
decaying extreme type II string. Our goal is calculating the barrier (potential energy 
versus 9) which will be later used in the calculation of the decay rate. We will need 
the kinetic term for the field 9(t, z) too, but we begin with the potential term. 

The condition ( p.l| ) ensures that both, the adjoint scalar and quark fields vary 
in small cores with sizes of order l/m a and l/m q , respectively. In the domain 
m~ l < r <C m^ 1 the fields 4> and q already reach their boundary values (pe(oo) = 
and qe(oo) = v. Moreover, we again use Eq. ( |5.2j) for the function fg, the leading 
logarithmic contribution to the tension coming from the boundary value fe(0) = 1. 

Substituting this data in Eq. ( fj.lip we get 



T l \6) = 2tt ( In ^) v 2 cos 2 9 + 2tt ( In — sin 2 29 . (5.4) 



m 7 ) \ m 7 1 



This is our result for the barrier profile for the extreme type II string. The first 
term here comes from the third term in Eq. ( |4.11| ) while the second one comes from 



the fifth term in Eq. (|4.11|) . All other contributions to T contain no large logarithms. 



At 9 = the tension in (|5.4| ) coincides with that of the ANO string, see Eq. (|5T 
At non-zero 9 the V 2 term dominates — it produce a huge barrier, with the height 
of order of V 2 In (m a /m 7 ). At 9 = tt/2 the tension T ll (8) vanishes, which means 
that the string disappears. This is summarized in Fig. ^ presenting T ll (8) versus 9. 

5.3 Effective world-sheet theory 

In order to calculate the decay rate of the string we have to work out the effective 
theory for the unstable mode parametrized by 9 on the string world-sheet. The 
collective coordinate 9 becomes a field of a 2D sigma model on the string world- 
sheet. The tension ( |5.4| ) gives us the potential term in the action of this sigma 
model. To complete the problem we need to know the kinetic term. 

In order to obtain the kinetic term we apply the standard strategy: we assume 
that 9 adiabatically depends on the world-sheet coordinates a. (Throughout the 
paper we use here the static gauge for the string in which o\ = t and o% = £3.) 
The kinetic term for 9 in the 2D sigma model comes from those in the action ( |2.1| ). 
To calculate the kinetic term for 9 it is convenient to rotate our field configuration 
( P~T ) into a "singular" gauge performing the gauge transformation with the matrix 



U 1 , see Eq. (|4.10|) . In this gauge the boundary values of the fields at infinity do 
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Figure 3: The potential energy T(8). 



not depend on 9 and give no contribution to the kinetic energy term, for instance, 
qk{x) = e k q e (r), and 

= V| +U~ 1 (A<j ) )U, (5.5) 

(remember, A0 falls off at infinity). The gauge field falls off rapidly at large dis- 
tances, 

Ai{x) =-t(d l U- 1 )Uf e (r). (5.6) 

We pause here to discuss a subtle point in the calculation of the kinetic term. 
With time dependence switched off, for the static string, the non-vanishing compo- 
nents of the gauge potential are A,-, (j = 1,2). The components A n with n = 0, 3 



vanish, see Eq. (|4.1| ). However, as soon as we allow 9 to depend on the world-sheet 
coordinates, the components A n with n = 0, 3 must become non-zero. Indeed, let us 
first assume that A n = (we will immediately see that this is a wrong assumption) . 
Consider the contribution of the kinetic term of the gauge field F?L. It is clear that 



the only piece contributes to the kinetic term of 9, 

F nj = d n A j -d j A n -i[A n ,A j ], j = 1,2, n = 0,3. (5.7) 
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If the components A n vanished, then one would obtain 

F nj = d n A j = {d n 9){d 9 A j ), (5.8) 
which, being combined with Eq. (|5.5| ) for the gauge potential Aj, yields 

fe(r) ■ (5.9) 



d n Aj = 2 (d n 6) Sji 



sin 26 — + cos 26 ( — sin a H — - cos a 



2 V2 2 



With this formula the expression for F 2 - is not even gauge invariant with respect to 
the gauge transformations depending on 6, see Eq. (|5.8|) . 

The fact that we made a mistake by assuming A n = manifests itself in the 
singularity in Eq. (|5.9| ) at r = (note that /e(0) = 1, see Eq. (|4.3|) ). 

Thus, in the calculation of the kinetic term of 6 one cannot avoid switching on 
the components A n with n = 0, 3 which, naturally, must be proportional to d n 6. 
The following ansatz for A n goes through: 

A n = -2 (d n 6) cos a - y sin a) a e {r) , (5.10) 

where ag(r) is a new profile function and the angle a is defined in Fig. [I]. Generally 
speaking, A n is parametrized by three distinct profile function accounting for three 
generators of SU(2). However, as it turns out, the single structure presented in 
Eq. ( |5.10| ), leads to a fully self-consistent and complete ansatz, with the r — > 
singularity in F n j cancelled. We do not need two other structures. 

Substituting Eq. ( |5.10| ) in Eq. ( |5.7| ) and ignoring terms that are a priori non- 
singular at r = 0, we get 



F n j = 2 (d n 6) Eje -| 



sin 26 — + cos 26 ( — sin a H — ^ cos a 



(l-Ofl(r)) , (5.11) 



where we assume that the gauge profile function fg — 1 + • • • at r <C l/m 7 . The 
reason why we can keep only the most singular terms is as follows. Let us remind 
that in Sect. |5.2| we found the potential term for the 6 field in the logarithmic 
approximation (i.e. the approximation in which only those terms are kept which 
contain large logarithms of the mass ratio). Our task in this section is to find the 
kinetic term for the 6 field in the very same logarithmic approximation. 

In order to cancel the actual divergence of / d 2 xF^j at r — > we must impose 
the following boundary condition: 

a e (0) = l, (5.12) 

as well as 

a e (oo) = 0. (5.13) 

In fact, for our purposes — determination of the kinetic term with the logarithmic 
accuracy — it is sufficient to use the step function model for ag(r) similar to that 
in Eq. Q, 

1 at r < R , 
ag{r) = { (5.14) 
at r > R , 
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where we introduced a new parameter R, to be determined below, (it is assumed 
that R < l/m 7 ). 

With the profile function ag(r) presented in Eq. ( |5.14| ), the contribution of the 
gauge term to the kinetic energy of the field 9 takes the form 



d 2 x — F a F^ a 



Rm~ 



(5.15) 



This is not the end of the story, however, since, in addition, we have to take into 
account the kinetic energy coming from the scalar kinetic terms in Eq. ( |2.1| ) (the 
second and the third terms). It is obvious that the third term, associated with the 
quark field, is proportional to v 2 and can be neglected as compared to the second 
term — the contribution of the adjoint scalar which is proportional to V 2 . As a 
result, using Eq. Q5.10 ) to calculate D n <f), we get for the total kinetic energy 

1 



\ In 



Rm, 



+ const ■ (V R 2 ) 



(5.16) 



The parameter R can now be determined from the requirement that the coefficient in 
front of (d n 8) 2 be minimal. Minimizing this with respect to R we find the condition 



1 



g 2 R 



V R. 



In other words, R turns out to be small, of the order of the inverse mass of the W 
boson, 

R — , (5.17) 

m w 

in full accord with what our physical intuition demands. With this value of the 
parameter R, the first logarithmic term in Eq. ( |5.16 ) dominates over the second one, 
which can be thus ignored with the logarithmic accuracy. As a result, combining 
together the kinetic term ( |5.16|) with the potential term (|5.4j ), we arrive at the 
following action of the 2D sigma model: 



S. 



str II 



2tt / (To- 



rn 



m w 



{drfif 



InZ^L) v 2 cos 2 9- [ln^ ] — sin 2 2fl 



V 2 



(5.18) 



This is our final result for the effective theory of the unstable #-mode on the 
string world-sheet in the primitive ansatz. In what follows we will assume that 
the masses mw and m a are of the same order of magnitude, i.e. A ~ g 2 . In the 
logarithmic approximation we can then replace the logarithms in the first and the 
third terms in Eq. ( |5.18| ) by a single logarithm, say lnmiy/m 7 . This assumption is 
by no means crucial. 
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5.4 Decays through tunneling ("bubbles of true vacua") 

To calculate the decay rate of the metastable string we may forget for a while about 
the microscopic theory fl2.1|) and turn our attention to the effective theory (|5.18|) . In 



this latter theory the string state is nothing but "the false vacuum state" at 6 f» 0, 
see Fig. |3], while the "no-string state" is the true vacuum at 6 ~ n/2. 

The metastable string decay occurs through the creation of the monopole-like 
objects: at a certain z = —z a magnetic charge is produced, accompanied by 
the production of an anticharge at z = zq, through tunneling. In the interval 
— zo < z < z the magnetic flux tube is eliminated. 

The corresponding process in the effective theory looks as follows. In the initial 
moment of time the theory resides in the false vacuum. Then it tunnels into the 
true one. The tunneling creates an interval of true vacuum, which subsequently 
experiences an unlimited classical expansion. 

The quantitative description of the false vacuum decay is well-developed within 
the quasiclassical approximation |23j, [25| , which is fully applicable to the 2D sigma 
model ( |5.18|) . The applicability of the quasiclassical approximation will become clear 
shortly. The parameter which regulates this approximation is V/v^> 1. 

Let us briefly review the general procedure [24, 25] of calculating the probability 



of the false vacuum decay (for a comprehensive review see Details (as well as 

the final answer) slightly depend on the space-time dimension. Since our effective 
world-sheet theory (|5.18f) is (1+1) -dimensional, we will focus on this case. 

An appropriate description of the tunneling probability implies a Euclidean ro- 
tation, 

t-^it. 

After the Euclidean rotation, the Euclidean action of the effective world-sheet theory 
takes the form 



m w 




S Eud = 2tt In ___ 1 / dV J " { Qf + {dz9 f + _ sin^ 26 



V 1 



m 7 J J ^ <T L ' ' J 2 

+ kv 2 cos 2 9] , (5.19) 

k = fln^Uln^V 1 . (5.20) 

In this formulation the false vacuum decay goes through the creation of a bubble of 
the true vacuum inside the false one, see Fig. f|. The bubble is a classical bounce 
solution in the potential V(8) = —T Il (8) which depends only on the radial variable 
\/t 2 + z 2 . In the thin wall approximation which is relevant to our problem the 
concrete form of the bounce solution in the model (|5.19 ) is not important. The only 



parameter we will need to know is the tension r of the bubble surface. 

The bounce solution has a negative mode associated with instability in the bubble 
size p. Integration over this negative mode produces an imaginary part of the 
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Figure 4: False vacuum decay through a bubble of true vacuum (in the Euclidean 
time) . 



vacuum energy. The latter determines the false vacuum decay rate 0, |25| which 
thus is proportional to 

r ~ exp(-S , bub bie), (5.21) 

where Stubbie is the classical action of the bounce. 

In the problem at hand the ratio of the critical size of the bubble to the bubble 
surface thickness is very large (it is regulated by the same parameter V/v ^> 1). 
Therefore, to calculate the bubble surface tension r one can neglect its curvature 
and consider a flat wall separating two vacua — one at 9 « and another at 9 m ir/2. 
Simultaneously we can (and should) neglect the term k v 2 cos 2 9 in Eq. ( |5.19| ) which 
is responsible for the non-degeneracy of these two vacua. With the term k v 2 cos 2 9 
switched off, the vacua at 9 ~ and 9 ~ it/2 become degenerate, and the flat wall 
perfectly stable. 

The tension of the flat wall is obtained from minimization of the energy functional 



r = 4vr(ln^ Jdz^ (d z 9) 2 + ^ sin 2 29 \ , (5.22) 
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with the boundary conditions 

7T 

9 = at z — > — oo , 9 = — at z — > oo . (5.23) 
The solution of the minimization condition is a well-known sine-Gordon soliton, 

9{z) = — jarcsin (tanhmiy^) + — | . (5.24) 
The mass of the sine-Gordon kink gives the wall tension r, 

r = 4vr(W) V -. (5.25) 

The wall thickness is inversely proportional to the W boson mass, 

A ~ ■ (5.26) 

The potential term in Eq. ( 5.19Q presents a huge barrier under which the false 
vacuum must tunnel. At the same time, the energy density difference between the 
false and true vacua — we denote it by £ — is small since it is determined by the 
third term in Eq. (|5.19 ), 

£^T^ = 2nv 2 \n (^j . (5.27) 



The action of the bubble on the tunneling trajectory is given by fl24 

SWbbic = 2vrp r - np 2 £ . (5.28) 
The size of the critical bubble p* is determined by the extremum of the action, 

p. = T - = 2(-) 2 (5.29) 
£ \ v J gv k 

The ratio of the bubble radius to the wall thickness is indeed very large, 

|~Q 2 »1, (5.30) 

which justifies the thin wall approximation. 

Substituting the critical size from Eq. (|5.29 ) in Eq. ( |5.28j ) we find the tunneling 
action, 5* = tit 2 /£. The probability of the false vacuum decay which is equal to 
the probability (per unit time per unit length) of the ANO string breaking is 



r bLaking ~ e S * = exp I — 7T — 



expj-^ 4 (m^H, (5.31) 
[ g z kv z \ m 7 J J 
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where the parameter k is defined in Eq. ( |5.20| ). We conclude this section by rewriting 
Eq. ( |5.31|) in terms of the ANO string tension ( |5.3|) , 



u . 16tt 3 / m w \ 2 V 2 . 

Tbreaking ~ eX P ^ - -— hi -— ) . (5.32) 

9 \ m 1 J 1 ANO 

We will see in Sect. ^|that the decay rate of extreme type I string, being expressed 
in terms of the type I string tension, is given by the very same formula. 

6 Extreme type I string 

In this section we will consider strings in the limit of very small quark masses. We 
still assume that the adjoint scalar mass is much larger than all other scales, in other 
words we impose the condition 

m w ~ m a > m 7 > m q . (6.1) 

We will begin with a brief review of the stable ANO strings in this case, and then 
turn to decaying strings. To avoid bulky expressions in what follows it will be 
convenient to introduce a parameter L, 

L = ln^, (6.2) 

m q 

and a function L(9), 

L{9) = In . (6.3) 

m q cos 

6.1 Type I ANO string 

Let us outline the solution (jTjJ for the ANO string in the Abelian Higgs model ( j2.7p 
under the condition m 7 3> m q . 

To the leading order in L = In (m 7 /m q ) the vortex solution has the following 
structure. The electromagnetic field is confined to a core with the radius R g which 
we will estimate momentarily. The profile function for the gauge field is given, 
approximately, by an expression similar to Eq. 



1 - 4 , r < R a 



f={ R °' (6.4) 
0, r > R g . 

Moreover, the quark field is close to zero inside this core. On the other hand, outside 
the core, the electromagnetic field is vanishingly small. At intermediate distances 

R g < r < — (6.5) 
m q 
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the scalar field satisfies the free equation of motion, see the first equation in ( |3.2j ), 
where the third and the last terms can be ignored. Its solution is as follows: 



<l(r) =v{l-(ln — ) (in ) } . ((3.(3) 

\ rm q J \ R g m q J 

At large distances, r ^> l/m q , the function q(r) approaches its VEV v, the rate of 
approach is exponential, 

q — > v (1 + C exp(— m q r)) . 



Now, substituting Eqs. fl6.4T ) and Q5.6|) in Eq. (|3.4j ) we get the tension of the 
static type I string function of _R„, 



^ANO — 27T 



const 2 / 1 



(6.7) 



The first term here comes from the first term in Eq. ( |3.4Q ) which is concentrated 
inside the core. The second term in Eq. ( |6.7|) comes from the logarithmically large 
region (|6.5|), where the quark field is given by Eq. (|6.6|). 

Minimizing the right-hand side of Eq. (|6.7|) with respect to R g we determine R g , 



1 1 / \ 2 

1 t0 I I m~ 



R g ~ — L* = — ln^I . (6. 

mi mi \ m n I 



This expression demonstrates that in the case at hand the size of the string core is 
logarithmically larger than l/m 7 , due to the presence of the light scalar q. 

Using this result for R g it is straightforward to evaluate the tension of the extreme 
type I string. To this end we plug Eq. ( |6. 8[ ) back in Eq. ( |6.7| ) and then obtain |15 



T^ NO = 2ttv 2 L-\ (6.9) 

where L is defined in Eq. ( p.2\ j. The tension T^ NO is saturated, in the logarithmic 
approximation, by the kinetic energy of the quark field (the second term in Eq. ( |3.4| ), 
the "surface" energy). All other terms in Eq. ( |3.4j ), as well as corrections to profile 
functions, yield contributions suppressed by extra powers of Lr x . 



6.2 Decaying type I string 

We now turn to the derivation of an effective 2D world-sheet theory for the unstable 
mode 8(t, z) for the extreme type I strings. Our first task is determination of the 
potential term. 

The description of the decaying type I string we are going to use runs parallel to 
that for type II, see Sects. |] and [5]. The adjoint scalar (fig varies from its boundary 
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value Q4.7|) to zero in a very narrow core whose size is of the order of l/m a . The 
profile function fg for the electromagnetic field is concentrated inside a larger core, 
of radius R g (9). Note that the parameter R g introduced above now becomes a 9 
dependent function. Inside this "electromagnetic" core the profile function fe is 
approximately given by Eq. ( |6.4| ) with R g replaced by R g (9). The quark field is 
very small inside this core, while outside it is given by Eq. (|6.6| ), again with the 
replacement R g — > R g {6). Assembling all these elements together and substituting 
in Eq. ( |4.11| ) we get 



T\6) = 2tt 



const — -— + v 2 In + In — — sm 2 29 

g 2 (R g (9)) 2 \ R g mJ \ mj 2 



(6.10) 

Next, for each given 9 one determines R g (9) by minimizing T l {9) with respect to 
R g . In this way one finds 

RM ~ —L, (6.11) 



(cf. Eq. ( |6.8| )). With this expression for R g {9) the tension of the decaying type I 
string versus 9 (in the logarithmic approximation) takes the form 



where the function L(9) is defined in Eq. ( |6.3|) . The boundary values are as follows. 
At 9 = the potential term T l {9) is equal to the static string tension (|6.9f) . At 
larger 9 the potential T l {9) develops a very high barrier (at the maximum T 1 (9) ~ 
V 2 ln(m a /m 7 )), and then it vanishes at 9 — tt/2. Qualitatively, the behavior of 
T l {9) is perfectly the same as that depicted in Fig. [| Equation ( |6.12| ) concludes 
our calculation of the potential term in the effective world-sheet action for the 9- 
mode of the type I string. 

Calculation of the kinetic term for the type I string repeats the same steps we 
made in Sect. |5.3| for type II and thus gives the same result for the kinetic term of 



type I as in Eq. ( |5.18|) . It is easy to understand why: the solution for the gauge 



field is essentially the same for the two cases. Assembling together the kinetic and 
potential terms we finally arrive at the following effective world-sheet theory for the 
unstable mode of the extreme type I string: 

s — 2 * / *° {7 ( h1 ^) (a - 9)2 - m - ( ln S t- 22 "} ■ < 613 ' 

We are now ready to consider the false vacuum decay in this sigma model. As 
in Sect. |5.4j, the decay rate is given by the formula 



r 2 " 



expl-7r-j, (6.14) 
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where £ is the difference between the energy densities in the false and true vacua, 
£ = ?a NO , see Eq. ( |6.9| ), while r is the tension of the flat wall separating the 
two vacua. Remember, to calculate the latter we neglect the term v 2 /L(8) in the 
potential energy. Hence, the wall solution, as well as r, are exactly the same as in 
Sect. |5T^ , see Eq. ( |5.25| ), provided that m a ~ mw, so that the logarithms of mw /m^ 
and m a /m 1 are the same. We accept this simplifying assumption. 

Substituting r and £ in Eq. (|6.14j ) we finally arrive at the decay rate (per unit 



length) of the extreme type I string, 

~ exp \-T ( ln ^) ( ln ^) } • (6 - 15) 

Needless to say that if we express this decay rate in terms of the type I string tension 
we get the same formula ( |5.32|) as for the type II string, 



r bre aking ~ exp ^ - — In — -j— \ . (6.16) 

[ 9 \ m 7 / 1 ANO J 

7 Strings in the theory with two adjoint matter 
fields 

In this section we consider the decay problem for strings in the theory with two 
adjoint fields, see Eq. ( |2.12j) . The light fields x surviving in the "macroscopic" 
QED limit, see Eq. ( |2.14| ), are written in the matrix form as follows: 



x = 4, , u + \ ■ (7.i; 




The ANO solution for the Z 2 string ( p.6|) corresponds to the minimal winding of the 
scalar matrix, 

X = -j= exp (-i y a) (® J J exp (i ^-a^j h 1/2 (r) . (7.2) 

In Fig. this trajectory corresponds to a semi-circle running along the equator and 
connecting two points of intersection of the equator with the 1-st axis. As we have 
already mentioned this string is stable — one cannot unwind it. 

However, strings with multiple winding numbers are metastable. Let us use the 
method developed in the previous sections to calculate the decay rate of the string 
with the winding number n = 2. The winding of the scalar field x i n Eq. ( |3.8| ) in 
the matrix notation takes the form 




X = -y= exp (-i r 3 a) [ i n ) exp (i r 3 a) h(r) . (7.3) 
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Here the scalar profile function satisfies the boundary conditions 

h(Q) = 0, h(oo) 



v. 



(7.4) 



To unwind this string we can use the same "unwinding" matrix ( |4.2j) which we 
used in the theory with the fundamental scalar. The "unwinding" ansatz takes the 
form 

X (x) = -^UnU^heir), 



A 



Aj(x) 







A 3 = 0, 



iUdjU-^l- f e (r)], j = 1,2 
VU^-U~ 1 + Acj), 



(7.5) 



where is given by ( |4.5|) . Substituting this ansatz in the action ( p.!2|) we get 
T{9) = 2tt J rdr j $L cos 2 9 + h' 2 + J± h 2 e cos 2 9 (l - i sin 2 



1 /2 1 
+ 2^ + 2^ 



<^(cos 29 - 2f g cos 2 5) + Vf e sin 20 



n 2 



V 2 ) 2 + a 



? 2 - 2V^e sin 20] 2 + jtffthl cos 2 20 }> . 



(7.6) 



Comparing this result with Eq. (|4.11 ), we see that the only terms which are modified 
are the ones associated with the light scalar \. In particular, the terms coming from 
the heavy scalar <ft , responsible for the huge potential barrier, stay intact. 

To proceed, let us consider the case of the extreme type II string assuming that 
the light scalar mass m x is much larger than the photon mass, m x ^> m 7 . Using 
essentially the same step-function model for the profile functions in ( [7.5|) as in Sects. 
|5.1| and [57|, we finally arrive at the following effective sigma model for the unstable 
mode on the string world sheet: 



S, 



str II 



2tt 



m w 



-4 In^i 



v 2 cos 2 9 1 



(djf 



sin 2 9 - In 



sin 2 29 



(7.7) 



At 9 = the potential in this sigma model reduces to the tension of n = 2 ANO 
string, 

(7- 



T-iII, n=2 o 2 a i m 9 

-Ta.no = 4 In — - 
m 1 
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The extra factor 4 here, as compared to Eq. (|5.3|), is due to the n 2 dependence of 
the extreme type II string tension on the winding number n. 

Now, to calculate the decay rate of this string we follow the same steps as in 
Sect. |5.4j . Note that the tension of the domain wall is still given by Eq. ( |5.22 ) and 



the only modification is due to the expression for E which is now given by the n = 2 
ANO string tension ( |7.8|) , 8 = T^^~ 2 . It is clear that, being expressed in terms of 
the string tension, the decay rate is given by the same expression ( |5.32|) , 



r breaking ~expj-^ln^j (7.9) 

We can calculate, with the same ease, the decay rate of the extreme type I string 
in the theory with two adjoint scalars starting from Eq. ( |7.6| ) and following the same 
procedure as in Sect. |6|. Clearly, the result for the decay rate, when expressed in 
terms of the ANO string tension T^q 2 , is given by the same formula ( |7.9|) . 

8 Lessons 

We pause here to summarize what we have learned from the calculation above and 
to abstract general features. 

The kink ( |5.24j ) represents a lump of energy, a bulge at the end of the broken 



string associated with the production of one unit of the magnetic charge, see Fig. ||. 



Magnetic 
charge 




m 



Unperturbed string 



m 



W 



Figure 5: The right-hand half of the broken string (in the ansatz ( |4.1|) ). One unit 
of the magnetic charge is produced in the shaded area. Arrows mark the magnetic 
field flux. 

The kink mass r (see Eq. ( |5.25| )) is the mass of the bulge carrying one unit 
of the magnetic charge. Keeping in mind that Mm ~ V/g one may ask why r is 
logarithmically larger than Mm- 

The answer to this question is quite obvious. While the ansatz (fO]) does describe 
the production of the magnetic charge at the end of the broken string, it is a highly 
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excited monopole-like state rather than the 't Hooft-Polyakov monopole. Indeed, 
the longitudinal dimension of the bulge is ~ rn^, a typical size of the monopole 
core. At the same time its transverse dimension (in the plane perpendicular to the 
string action) is of order m" 1 . This is much larger than the monopole core size. The 
stretching of the core in the perpendicular direction is the reason why this lump is 
logarithmically heavier than the 't Hooft-Polyakov monopole. This is an inevitable 
consequences of the fact that the ansatz ( |4. 1|) contains a single profile function fe{v) 
which governs the behavior of both, the photon and the W boson fields. 

By introducing two distinct profile functions one can readily eliminate the short- 
coming in the ansatz ( f4.1| ). We do it in Appendix. However, even before particular 



improvements are done we want to note that sufficient experience is already ac- 
cumulated to enable us to develop a proper effective description of the tunneling 
process. 

Indeed, in actuality the end-point domain of the broken string is roughly a hemi- 
sphere with the radius ~ vn~ x . The core which emanates the magnetic flux has 
dimension of order of m^ 1 in both transverse and longitudinal directions. In fact, 
the core is practically unperturbed 't Hooft-Polyakov monopole. This is due to the 
fact that at distances of order the effect of the (magnetic charge) confinement 
is negligible, it comes into play only at distances ~ m~ l . Thus, the mass of the 
end-point bulge in a "good" ansatz must be r = M M + 0{y/g). The 0{y/g) correc- 
tion reflects the distortion of the 't Hooft-Polyakov monopole at distances ~ m' 1 , 
and can be neglected compared to Mm- A potential of the type depicted in Fig. |3| 
will emerge leading to the tunneling problem described by the bubble action (see 
Appendix for details), 

^bubble = 2tt pM M - np 2 T ANO . (8.1) 



The only distinction with the consideration carried out in Sect. |5.4j is that in the 
"good" ansatz the ratio of the bubble radius to its thickness 



P* (V 
A \v 



(8.2) 



(cf. Eq. ( |5.3U| )). We loose one power of V/v, but the ratio p*/A is still large, and 
the thin wall approximation justified. 

The action ( |8.1| ) immediately leads to the string decay rate given in Eq. ( |1.3| ). 
The effective approach based on Eq. ( |3.1|) is exactly the one of Vilenkin and 
Preskill and Vilenkin 



In view of a general nature of this conclusion, let us comment on the hierarchy 
of parameters we deal with. We want to separate elements of this hierarchy which 
are absolutely essential for our consideration from those which bear a technical 
character. 

Our consideration is quasiclassical. It is valid only at weak coupling. This 



requires Eq. (1.1) to be valid. We cannot sacrifice this condition. 
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Moreover, the two-stage nature of the gauge symmetry breaking, Eq. ( |1.2|) , is 
crucial too. Among other things, it guarantees that the thin wall approximation is 
justified. It is unclear whether lifting the condition ( |1.2[ ) one can still develop an 
analytic description of the tunneling process. We will not try to lift the constraint 

On the other hand, the relation between m q and m 7 is clearly of a technical 
nature. In the two extreme cases, when m g /m 7 is either very large or very small 
one can calculate ANO string tension analytically. One can pose a question what 
happens at arbitrary values of m q /m T Calculation of the ANO string tension in this 
case certainly requires numerical computations. However, this calculation can be 
carried out entirely in the macroscopic low-energy U(l) theory, with no reference to 
the microscopic non-Abelian theory. The procedure is well-developed; the function 
F(m g /m 7 ) which parametrizes the ANO string tension in the general case, 



Tano = 2t™ 2 F(^) , (8.3) 



is known in the literature 32 . 



Irrespective of the value of m g /m 7 , the energy difference in the false and true 
vacua is £ = Tano Moreover, the expression ( |1.3|) for the decay rate is universally 
valid as long as p*/A ^> 1, i.e. V/v ^> 1. The flat wall tension r is calculated in the 
approximation which neglects the quark field contributions altogether — the barrier 
is determined only by terms proportional to V 2 . No matter what the particular 
ansatz is, it must yield r = Mm modulo possible small correction ~ v/g. This 
condition can be viewed as a test of the "goodness" of the ansatz. 

Therefore, if one neglects all terms suppressed by powers of v/V one inevitably 
arrives at Eq. (|3.1|) , irrespective of the ratio m q /m T The result ( |L3"1) for for the 
string decay rate r breaking ensues. It is valid for arbitrary value of the ratio m g /m 7 . 

Of particular interest are examples emerging in the supersymmetric setting. For 
instance, for the BPS string Tano — 27rt> 2 . Moreover, the BPS bound for the 
monopole mass Mm is 

M M = ^^V. (8.4) 
9 

Then Eq. (|1.3j ) predicts the following decay rate of this string: 

f 8n 2 2V 2 ) , . 

r BP s~exp|-— — | . (8.5) 

Let us recall that the Abelian BPS strings embedded in non-Abelian gauge theories 



appear, say, in the charge vacua of the Seiberg-Witten theory with matter |TJL 129 



30, 17 . 
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9 Comparison with Schwinger's expression 



The decay of the string goes through breaking of the string in pieces through a 
monopole-antimonopole pair production. From this standpoint, the bubbles of the 
"true vacuum" of the effective 2D sigma model inside the "false" one are, in fact, 
domains where the string is broken by monopole-antimonopole pairs. The decay 
rate is exponentially small for large monopole masses; the exponent is determined 
by the ratio M^/Tano 

Given this interpretation it is instructive to compare Eq. (|1.3|) with the famous 



Schwinger's formula |2(J, |2T| for the electron-positron pair production in the constant 
electric field, 

where m is the electron mass, and E is the electric field 0. Note that with our 
normalization, see Eq. Q2.1|) , the coupling constant is included in E. The probability 
fl9.1| ) can be obtained as the imaginary part of the one-loop graph presented in Fig. |6|. 



Figure 6: Propagation of the electric charge in the background electric field. The 
charged particle loop is denoted by a thick line, the background field by dashed 
lines. 

Dualizing Schwinger's expression we can try to use it for evaluating the proba- 
bility of the monopole-antimonopole pair creation in the magnetic field existing in 
the core of the ANO string. Of course, we have to assume that this field is constant 
on the scale of the monopole-antimonopole separation (we hasten to add that this 
is a wrong assumption). 

It is not difficult to get from Eq. ( |9.1| ) a dualized Schwinger formula for the 
magnetic monopole pair production in the homogeneous magnetic field. Indeed, 
with our normalization the duality transformation reads E <-» g~ 2 B where B is the 

2 Extensions of the Schwinger formula in string theory were recently discussed in Refs. pi B4|. 
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magnetic field. Then the dualized Schwinger formula takes the form 



2 



Tmm ~ exp (-„ . (9.2) 

The magnetic field in the ANO string can be readily estimated from its flux, 

2tt 1 

B m 2 = -g 2 v 2 . (9.3) 

area 1 2 



Combining Eqs. (|9.2| ) and (9.3) one obtains 



r M M~exp I-ttC-^ , (9.4) 



v 2 



where C is a numerical coefficient of order 1 . 

One might ask whether dualizing Schwinger's formula is a good idea. Indeed, 
Schwinger 's formula assumes the smallness of the gauge coupling constant g 2 . The 
coupling of the magnetic monopoles ~ g~ 2 is then strong. Therefore, the 7 quanta 
exchanges, as in Fig. |7|, might drastically change the result. 

This effect was studied in Ref. |21|| , where it was found that the extra term in 



the exponent due to the gamma quanta exchanges is of order 1/g 2 . This is easy to 
understand. Indeed, the corresponding Coulomb energy is of order l/(p*g 2 ) which 
generates a correction in the bubble action of the order of 1/g 2 . It can be safely 
neglected because M M /v 2 3> 1/g 2 . 

Comparing this Schwinger- formula-based expectation with our result ( |1.3[ ), and 
keeping in mind that Tano — 27rt> 2 modulo a logarithm of m q /m 7 , we see that 
the powers of the mass scales and of the coupling g 2 in the two exponents per- 
fectly match each other. However, the logarithmic factors (in are missing 
in ( |9.4| ). These will not (and need not) match, and neither a numerical factor in 
front of M M /v 2 . The reason for this is that Schwinger's formula and our calcula- 
tion refers to physically different phases of the theory. Schwinger's formula ( |9.2|) 
describes monopole-antimonopole production in the external magnetic field in the 
Coulomb phase while our result ( p..3|) refers to the breaking of the ANO string by 
the monopole-antimonopole production in the confinement (for monopoles) phase 
of the theory. 



10 Conclusions 

In this paper we calculated the decay rate of an Abelian flux tube embedded into a 
non-Abelian theory. We focussed on two simplest examples where the phenomenon 
does occur: non-supersymmetric SU(2) gauge theory with the adjoint and funda- 
mental scalars and the same theory with two adjoint scalars. In the first example 
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Figure 7: Propagation of the monopole in the background magnetic field. The 
monopole loop is denoted by a thick line. 7 quanta (denoted by a zigzagy line) are 
strongly coupled to the monopoles, and exchange of the 7 quanta must be taken 
into account. 



all ANO strings are metastable. In the second one the string with minimal wind- 
ing number is stable (Z 2 -string), while strings with multiple winding numbers are 
metastable. 

The pattern of the SU(2) symmetry breaking is two-stage, 

SU(2) U(l) nothing, V > v . 

As was expected, the decay rate is exponentially small. The suppressing exponent 
is proportional to the ratio of monopole mass squared to the string tension. The 
interpretation of this result is that the string gets broken into pieces by the monopole- 
anti-monopole production. Note that the monopole-anti-monopole pair production 
in the given context by no means implies that superconductivity is lost. 

Although we considered a particular model with metastable strings, we believe 
that the final answer is rather general and can be qualitatively applied to any 
metastable Abelian string embedded in a non-Abelian theory. In particular, as 
we mentioned in Sect. [I], the reduction of string multiplicity from Z N ~ X down to 
in the strong coupling vacua of the Seiberg-Witten theory is due to a similar mech- 
anism. In this case we deal with electric strings (which arise due to the monopole 
/dyon condensation), so that the metastable strings must broken by the W boson 
pair production (rather than the monopole pair production which takes place in the 
magnetic flux tubes). 

Unwinding ansdtze of the type presented in Eq. (|4.1D can be used in other similar 
problems, for instance, for studying the metastability of the appropriately embedded 
semilocal strings (cf. f23j, for a review of the semilocal strings see Ref. f35fl). 



Finally, it is worth noting that the calculation of the string decay rate presented 
here can be viewed as a calculation of an open string coupling constant in the 
effective string theory of ANO string. 
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Appendix. Improving the primitive unwinding 

ansatz 



As we explained in Sect. || our ansatz (|4.1|) gives a monopole in a highly exited 
state at the end of the broken string. Its mass r is much larger then the monopole 
mass due to the logarithmic factor in (|5.25|) . The reason is that we used the same 
profile function fg for both a = 3 and a = 1, 2 components of the gauge potential 
in ( [4.1| ). We used the model Q5.2| ) for this profile function which ensures that all 
components of the gauge field are nonvanishing inside the region of size ~ l/m 7 
in the (l,2)-plane. However, it is clear that the a = 1, 2 components of the gauge 
field are very heavy (W bosons) and should be spread over a much smaller region, 
of size ~ 1/mw Now, let us modify the ansatz (|4.1| ) to take this circumstance into 
account. 

First, let us explicitly write down the gauge field given by the ansatz ( |4.1|) in the 
singular gauge, see Eq. ( |4.10| ), 



Ai 



— e 



V r 2 



2 cos 2 29 



T 3 



sin 29 



n . r 2 

— sin a H cos a 

2 2 



fe(r) 



(A.l) 



To modify it, we introduce two different profile functions, in front of r 3 and r 12 
matrices, as follows: 



Ai 



x 



2cos 2 2# ^f e \ 



r) — sin 26 



— sin a H cos a 

2 2 



fi 



w , 



(A.2) 



We assume that / is concentrated inside a very small region of size ~ 1/mw, 
whereas / 7 is spread over a much larger domain of size ~ l/m 7 . The expressions 
for the adjoint and the fundamental scalars are still given by Eq. fl4.10|) . 
Substituting this into the action (|2.1|) we end up with 



T{B) 



cos 2 9 




cos 2 W)' 2 + sin 2 ^ 



W\i 2 



+ 



/ 2 
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+ 



Te 



cos 2 9 (cos 2 9f e + sin 2 Qff 



+ 



+ 



1 , 



+ 



2r 2 



\p e (cos 29 (1 - f e ) - (cos 2 2#/ e 7 + sin 2 2dff) ) + V/f sin 20] 2 



+ A 



2VVe sin 26* 



2 7 

+ 4 



-2 2 



(A.3) 



We see that now the term V fg sin 29 in the square brackets contains f w rather 
than / 7 (cf. Eq. ( |4.11| )). This means that its contribution in the integral will be 
saturated in a a small region, of size ~ 1/mw, and will not lead to large logarithms. 
Thus, the height of the barrier now is V 2 instead of V 2 ln(m a /m 7 ). It should be 
added, however, that in the absence of large logarithms we cannot obtain the barrier 
profile using analytic calculations. We need computer simulations to minimize the 
tension (|A.3|) and find the profile functions. This goes beyond the scope of this 
paper. 

Still, some simple examples of the profile functions we have analyzed show that 
we can avoid getting large logarithms both in the barrier height and in the kinetic 
term for the 9 field in the effective 2D sigma model for the unstable string mode. 
This leads us to the following representation for the domain wall tension r 



T 



const 



J dzjl (d z 9) 2 + ^-sm 2 29^ , 



(A.4) 



instead of the one in Eq. ( p. 22] ). Here, the constant in front of this action can be 
fixed by means of a numerical minimization procedure in ( |A.3j ). Note, that the 
relative coefficient between the kinetic and the potential terms in Eq. (|A.4 ) is fixed 
by the requirement that the mass of the 9 field should coincide with that of the W 
boson. 

Calculating the domain wall tension in the theory ( |A.4|) we get 



const 



V 
9 



(A.5) 



In order to prove that actually 



r Mm (A.6) 
we need numerical simulations. Still, if we accept on physical grounds that this is 



the case, then the bubble action in the 2D sigma model given by Eq. flS.ip ensues. 
This leads us to the final formula ( |1.3| ) for the decay rate of the metastable string. 
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